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Based on fundamental particle physics processes like the production and subsequent decay of 
pions in interactions of high-energy particles, close connections exist between the acceleration 
sites of high-energy cosmic rays and the emission of high-energy gamma rays and high-energy 
neutrinos. In most cases these connections provide both spatial and temporal correlations of 
the different emitted particles. The combination of the complementary information provided by 
these messengers allows to lift ambiguities in the interpretation of the data and enables novel and 
highly sensitive analyses. 

In this contribution the H.E.S.S. multi-messenger program is introduced and described. The 
current core of this newly installed program is the combination of high-energy neutrinos and 
high-energy gamma rays. The search for gamma-ray emission following gravitational wave 
triggers is also discussed. Furthermore, the existing program for following triggers in the 
electromagnetic regime was extended by the search for gamma-ray emission from Fast Radio 
Bursts (FRBs). An overview over current and planned analyses is given and recent results are 
presented. 
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1. Introduction 

The key question to resolve the long standing mystery of the origin of high-energy cosmic 
rays is to locate the sources and study the acceleration mechanisms able to produce fundamental 
particles with energies orders of magnitude above man-made accelerators. The combination of 
complementary information provided by multiple messengers and novel techniques will increase 
the chances to achieve this century old task. 

Thanks to the increase in computing power combined with advances in electronics and data analysis 
techniques, real-time analyses of data taken by astrophysical observatories at all wavelengths are 
now becoming possible. In addition, new messengers provide additional information on high- 
energy astrophysical processes. Due to fundamental particle physics processes like the production 
and subsequent decay of pions in interactions of high-energy particles, the acceleration sites of 
high-energy cosmic rays might also be sites of gamma ray and neutrino production. Unlike cosmic 
rays, both gamma rays and neutrinos are neutral messengers and therefore point back to their 
source. If produced within the accessible horizon (which is limited for gamma rays due to pair 
creation processes on the extragalactic background light) spatial and temporal correlations of the 
different emitted particles could exist. The H.E.S.S. multi-messenger program will fully exploit 
these correlations and might allow opening a new window to the high-energy universe: real-time 
multi-messenger astronomy at TeV energies. 

2. High-energy neutrinos 

The core of the H.E.S.S. multi-messenger program is to exploit the intimate connection be¬ 
tween high-energy neutrinos and y-rays. Provided appropriate conditions of the environment of 
cosmic accelerators (e.g. magnetic helds, matter and held densities, etc.), high-energy (hadronic) 
particles are potentially undergoing interactions with matter and radiations helds within and/or 
surrounding the acceleration sites. The light mesons, predominately pions, created in these inter¬ 
actions will decay by emitting both high-energy neutrino as well as y-rays. 



Figure 1: Arrival directions in Galactic coordinates of high-energy neutrino events detected by IceCube. 
Here we study the regions around events IC-5 and IC-18 (highlighted by the black circles) . Modihed 
from [4]. 
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Two major instruments searching for these astrophysical neutrinos are currently in operation: 
IceCube at the South Pole [1], complemented by ANTARES [2], a smaller detector in the Mediter¬ 
ranean Sea. Despite tremendous efforts, so far none of the neutrino telescopes has found any 
significant localized excess (e.g. [3]). Yet, a significant breakthrough in the search for astrophys¬ 
ical neutrinos has been made very recently by the IceCube collaboration. In three years of data, 
IceCube was able to single out 37 neutrinos with energies in the range of 30 TeV to 2 PeV that 
interacted within the instrumented volume (see Fig. 1). The atmospheric background contribution 
has been estimated as 8.4 zb 4.2 from cosmic ray muon events and 6.6^^'^ atmospheric neutrinos 
events, which yields a statistical significance of 5.7(7 [4]. However, the origin of these neutrinos is 
unknown and no significant clustering or excess at small angular scales has been found so far. 

Here we present detailed studies of the regions around some of the neutrino events detected 
by IceCube, focussing on track-like events that provide a good angular resolution of the order of 
1°. The studied regions are highlighted in Fig. 1. The main aim of the presented analysis is the 
search for steady gamma-ray emitters in the vicinity of the arrival direction of the observed high- 
energy neutrino events, the discovery of which would trigger detailed follow-up studies at various 
wavelengths and messengers to clarify a potential association. The two selected regions around 
events with the identifiers IC-5 and IC-18 (see [4]) have been studied in the GeV y-ray domain 
using data from the the FAT instrument onboard the Fermi satellite and data taken by the H.E.S.S. 
Imaging Air Cherenkov Telescopes (lACTs), which are sensitive to multi-GeV-TeV y-rays. 

2.1 H.E.S.S. analyses 

Dedicated observations on the two IceCube events IC-5 and IC-18 were taken with the H.E.S.S. 
lACTs during 2014. 

The region around IC-5 has been observed for about 2 h at a zenith angle of 30° in monoscopic 
mode with the H.E.S.S.-II 28m telescope. After correcting for acceptance effects the effective live 
time corresponds to 0.9 h. The data were analyzed using the Model Analysis [5] with standard 
gamma-hadron separation and event selection cuts and using a background estimation exploiting 
the uniformity of the acceptance across the field-of-view (FoV) of the system. No significantly 
enhanced gamma-ray flux has been detected. The distribution of gamma-ray events exceeding the 
background is shown for the full ROI in the left plot of Fig. 2. The right plot of Fig. 2 shows the 
map of the Li & Ma significances [6]. They are both fully compatible with the background expec¬ 
tation. 

The region around the neutrino event IC-18 has been observed for almost 12 h (effective live time 
after acceptance correction of 9.5 h) with the full array of the five H.E.S.S. telescopes at zenith 
angles around 10°. The data were again analyzed using the Model Analysis [5] with standard 
gamma-hadron separation and event selection cuts. Thanks to the high available statistics the back¬ 
ground could be determined using the "reflected background" method described in [7], a method 
that exploits the properties of the wobble data taking mode and yields very low systematic uncer¬ 
tainties related to the acceptance of the camera system. No significantly enhanced gamma-ray flux 
has been detected and the resulting maps (cf. Fig. 2) are fully compatible with the background 
expectation. 
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Figure 2: H.E.S.S. observations of IceCube events 5 and 18. The dashed circles denotes the individual 
median angular uncertainty on the neutrino directions provided by IceCube [4]. For comparison, the inner, 
solid circles show the approximate size of a point-like source after convolution with the H.E.S.S. angular 
resolution. Left plots: Map of VHE gamma-ray events exceeding the background expectation. Right plots: 
Point source significance map. 


2.2 Fermi-LAT analyses 

Fermi-LAT data recorded between 2008-08-04 and 2015-05-01 has been analyzed (P7Rep) in 
the full available energy range (100 MeV - 300 GeV). All photon candidates (evclass=2) above 
£* = 100 MeV fulfilling the basic quality criteria proposed by the Fermi collaboration (quality=l; 
LAT config=l; rock angle< 52°; zenith< 100°) within a 10° x 10° region of interest centered at 
the two locations of interest (IC-5: RA = 110.6°,Dec = —0.4° and IC-18: RA = 345.6°,Dec = 
—24.8°) have been selected. The count maps of the selected gamma-ray events is shown in the left 
column of Fig. 3. The gamma-ray emission observed in the Fermi-LAT data has been modeled 
using the information given in the 3FGL catalogue [8], fixing the parameters for known sources 
further away than 3 deg from the center of the region of interest (ROI) to the 3FGL values. An 
additional point-like source with a power-law energy spectrum has been added to the description 
before fitting the model parameters to the count map of the selected events. The fit did not yield 
significant emission from the additional putative source in either of the studied fields (cf. Fig. 3). It 
should be noted that the ROI of IceCube-5 (median angular resolution of 1.2°) comprises an AGN 
know to emit y-rays PKS 0723-008 which is located 1.04° from the center of the ROI. Given the 
high chance probability of finding a known y-ray emitting AGN within the error-box of a track-like 
high-energy neutrino detected by IceCube of about 37 % and based on flux level and spectrum of 
PKS 0723-008, an association with IceCube-5 can be discarded (e.g. [9]). 
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Figure 3: Region around the studied IceCube events as seen by the Fermi-LAT gamma-ray observatory. 
Left plots: Count maps of y-rays detected by Fermi-LAT in equatorial coordinates. Center plots: Maps 
showing the model of y-ray sources fitted to the Fermi-LAT data. Right plots: Normalized residual maps. 
The white circles denotes the individual median angular uncertainty on the neutrino directions provided by 
IceCube [4]. 


2.3 Upper limits on the y-ray flux 

Given the absence of a significant very high energy y-ray signal in the observed regions upper 
limits on the VHE and HE y-ray flux have been derived. Elux limits <I>ul have been calculated 
for a point-like source in the center of the ROI (see solid circles in Eig. 2) and for a generic E~^ 
energy spectrum. Eollowing the method introduced by Eeldman & Cousins [10], the obtained 99 % 
confidence level limits from the H.E.S.S. observations are shown as black triangles in Eig. 4. The 
limits on the high energy y-ray flux derived from the Eermi-LAT data are shown as arrows. Eor 
reference, the red markers denote the expectation from the neutrino candidate events. The esti¬ 
mation of the neutrino flux related to the individual events (dashed markers) follows the approach 
described in [11], assuming that each single event comes from one astrophysical counterpart with 
the observed flux following an E~^ energy spectrum and being spread over 1 dex in energy. The 
conversion into a y-ray flux (solid red markers) relies on Monte Carlo simulations of the hadronic 
interactions connecting neutrino and y-ray fluxes via the decay of charge and neutral pions cre¬ 
ated in pp interactions within or close to a generic hadronic accelerator. This conversion relies 
on several assumptions on the source environment like sufficiently low radiation and matter densi¬ 
ties to prevent y-ray absorption, pion-matter and ;?y interactions. Details of these assumptions and 
considerations are given in [12]. 

Unlike neutrinos, high energy y-ray photons are absorbed by pair production on the extra- 
galactic background light (EBL). This process can be described by ^obs = ^source x ^ where 
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Figure 4: VHE gamma-ray flux limits ^ul at 99 % CL derived from the H.E.S.S. (black triangles) and 
Fermi-LAT (black arrows) observations assuming a point-like source with an E~'^ energy spectrum. The 
estimate of the y-ray flux (solid, red marker) derived from the IceCube observation (dashed red marker) of 
event IC-5 (left plot) and IC-18 (right plot) is shown for reference. 

the optical depth T is a function of the energy Ey and the redshift of the source Zs- For sufficiently 
distant sources, i.e. sufficiently large optical depths, the expected gamma-ray flux <t>y will get ab¬ 
sorbed and could therefore become compatible with the upper limits <I>ul derived from the H.E.S.S. 
measurements. Using the EBL model given in [13], we derive these minimal distances to z = 0.007 
(0.012) for IC-5 and IC-18 respectively. It should be noted that this calculation very conservatively 
assumes that the potential sources are only emitting in the 100 TeV range around the energies mea¬ 
sured by IceCube, i.e. no lower energy radiation is present. Extrapolation of the IceCube flux to 
lower energies would allow to put severe constrains on the source distance, but would induce a 
significant dependence on the assumed spectral shapes. 

2.4 Real-time neutrino alerts 

Despite several efforts, including the ones presented here, the origin and astrophysical coun¬ 
terparts of the high-energy neutrinos detected by IceCube remain unknown. This hints to the pos¬ 
sibility that the underlying sources are either very faint but numerous or of transient nature and 
emit neutrinos only for a limited amount of time. Should the latter be the case, the time domain 
of multi-messenger and multi-wavelengths searches becomes increasingly important. Only by trig¬ 
gering deep follow-up observations of significant neutrino events rapidly after their occurrence one 
can be sure to obtain complete contemporaneous multi-messenger and multi-wavelength coverage, 
which might be necessary for the unequivocal proof of a common origin of potentially observed 
transient events. 

The H.E.S.S. multi-messenger program is actively pursuing this direction. Over the last years 
the second phase of the H.E.S.S. experiment has been commissioned successfully. It now includes 
the H.E.S.S.-II 28m telescope, the largest Cherenkov telescope in the world and providing the low¬ 
est energy threshold of ground based gamma-ray detectors worldwide. Another major goal for 
this new phase of H.E.S.S. is to reduce the response time of the system in order to increase the 
capabilities for the detection of transient phenomena. As discussed in [14], this goal has been 
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Figure 5: Localization probability of the origin of a simulated gravitational wave detected jointly by the 
advanced Virgo and Ligo interferometers (from 1. Bartos [16]). 


achieved on the level of the instrument conception and construction. Significant further improve¬ 
ments to the DAQ and alert system have been introduced recently [15]. The search for transient 
multi-messenger signals will greatly benefit from these developments. For example, close collab¬ 
orations between H.E.S.S. and the two major neutrino telescopes have been established and will 
enable rapid follow-up observations of high-energy neutrino candidates to start in the very near 
future. 


3. Gravitational waves 

In addition to the above-mentioned recent discoveries of extra-terrestrial, high-energy neu¬ 
trinos, a completely new messenger might soon become available: gravitational waves. The two 
major, kilometer-scale interferometers searching for gravitational waves, Ligo and Virgo, are about 
to finish extensive upgrade programs and will start taking physics data in the next years (Ligo in 
2015 and Virgo in 2016/2017). Predicted by the theory of General Relativity, direct detection of 
gravitational waves and the localization of their origin will permit important studies of both the 
underlying emission mechanisms as well as tests of fundamental principles of physics. The main 
signals for the advanced detectors are expected to come from the coalescence of binary neutron 
star systems. As these systems are also thought to be the progenitors for short gamma ray bursts, 
a connection between gravitational waves and other messengers like neutrinos, cosmic rays and 
gamma rays might be possible. 

Unfortunately the gravitational wave detectors will not be able to precisely localize the origin 
of the detected emission. Fig. 5 shows a typical example for the localization probability of the 
origin of a gravitational wave detected jointly by Virgo and Ligo. An extensive follow-up program 
involving a large range of instruments is therefore being set up. Thanks to their high sensitivity 
and large FoV, high-energy gamma ray observatories will play an important role in this endeavor 
(e.g. [17]). The transient nature of the expected events (e.g. afterglows of short GRBs fade away 
within minutes) makes rapid follow-up responses a crucial requirement. Thanks to its rapid re¬ 
sponse time (combined with its high sensitivity and large FoV), the H.E.S.S. y-ray observatory is 
well suited for these observations and preparations for the first physics runs of the advanced GW 
detectors are in progress. 
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4. Multi-frequency transients 

H.E.S.S. continues to pursue its long-standing follow-up program on triggers from observato¬ 
ries covering the full electromagnetic spectrum from radio to high-energy y-rays. For example, a 
dedicated program is performing follow-up observations on Gamma-Ray Bursts [15]. 

This program has been extended recently to include radio transients and especially Fast Radio 
Bursts (FRBs). Several of these extremely powerful and rapid bursts have been discovered in 
archival radio data over the last past years [18]. They are considered to be emitted by extra-galactic 
objects but the underlying process remains elusive. Due to their brevity (lasting only a few mil¬ 
liseconds), online detections and rapid alerts to follow-up facilities play again a crucial role in the 
quest for their origin. Several radio observatories are currently putting in place such data analysis 
pipelines and first follow-up campaigns are being performed. H.E.S.S. is currently collaborating 
with the SUPERB program at the Parkes radio observatory, able to emit FRB alerts with a delay 
O(20s) [19]. 
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